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It is shown how it is possible to reconstruct the initial state of a one- dimensional system by 
measuring sequentially two conjugate variables. The procedure relies on the quasi-characteristic 
function, the Fourier-transform of the Wigner quasi-probability. The proper characteristic function 
obtained by Fourier-transforming the experimentally accessible joint probability of observing "po- 
sition" then "momentum" (or vice versa) can be expressed as a product of the quasi-characteristic 
function of the two detectors and that, unknown, of the quantum system. This allows state re- 
construction through the sequence: data collection, Fourier-transform, algebraic operation, inverse 
Fourier-transform. The strength of the measurement should be intermediate for the procedure to 
work. 



Introduction. Quantum state tomography, i.e., the re- 
construction of the unknown state of a quantum system 
is a fundamental problem. Its formulation can be traced 
back to Pauli [1], who asked whether a measurement of 
position on an ensemble of spinless systems prepared in a 
pure state, complemented by a measurement of momen- 
tum on a distinct ensemble, would allow to reconstruct 
the wave function. This question has been answered in 
the negative [2, 3]. Pure states are represented by a wave- 
function, and they are exceptional, in the sense that they 
form a zero-measure subset of all possible states (which 
reflects in the great care one has to take in preparing a 
pure state). The most general state is mixed and usually 
described by a density matrix, which is then the object 
to be reconstructed. However, other equivalent descrip- 
tions of a mixed state have revealed to be more useful or 
significant than the density matrix: the Wigner function 
[4], the Husimi Q function [5], the Glauber-Sudarshan 
P function [6, 7], the recently reintroduced Dirac func- 
tion [8], and the almost forgotten Fourier-transform of 
Wigner function, which we refer to as the Moyal M func- 
tion [9]. Furthermore, parametric families including all 
of the above functions have also been introduced [10, 11]. 

Despite the quantum state being essential in describ- 
ing a system, successful quantum state reconstruction, 
dubbed quantum state tomography because of an anal- 
ogy with the germane procedure of medical tomography, 
is relatively recent [12]. The procedure relies on a pro- 
posal by Vogcl and Risken [13]. For a recent review 
of continuous-variable quantum state tomography, the 
reader may refer to Ref. [14]. 

Recently, a remarkable experiment [15] showed how 
it is possible to determine the unknown pure state of 
a one-dimensional quantum system by making a weak 
measurement of the x variable followed by a strong mea- 
surement of the conjugate variable p. A method, based 
on the Dirac function, allowing to lift the restriction to 
pure states was proposed recently [16]. Here, we make 
an alternative proposal, allowing the reconstruction of 
the Moyal M function by a quick sequence of two mea- 



surements of conjugate variables. In order for the pro- 
cedure to work, the strength of the measurement is not 
a fundamental issue, provided it is not too strong nor 
too weak. The equations at the basis of this proposal 
were reported, without derivation, in Ref. [17], and they 
are an exact result, not a perturbative expansion in the 
couplings. Furthermore, the procedure proposed herein 
requires a fixed setup. This is to be contrasted to ho- 
modyne quantum state tomography, where a different 
quadrature operator x$ is measured for various values 
of <fi. The efficiency of the two procedures is otherwise 
comparable, as in the sequential measurements one has 
to evaluate a joint probability of two variables, while in 
the homodyne detection scheme a one-parameter family 
of single variable probabilities is estimated. 

Review of the measurement. A linear measurement re- 
lies on a specific bilinear coupling between a quantum 
system and a probe. The interaction is assumed bilinear 
in the observable X belonging to the system and in an 
observable $ belonging to the probe: 



(1) 



where the function g(t) is strongly peaked around t = 
and has unit integral. If the spectrum of X is bounded, 
then A may include a scale such that the eigenvalues of X, 
indicated by X, are dimensionless and less than one [18]. 
A detector with a continuous unbounded output is con- 
sidered, so that $ has a conjugate operator J, satisfying 
[$, J] = i. Notice that XX and J are homogeneous. The 
variable J represents eigenvalues of J and is the readout 
of the detector, carrying information about the system. If 
the probe is initially in a well defined state with vanishing 
variance in the readout variable J, i.e. its density ma- 
trix is p(J, J') = 6j < j'6(J), then the measurement is an 
ideal strong one (we are indicating by S a ^ the Kronecker 
delta and by 6(a) the Dirac delta). When this latter re- 
quirement is relaxed the measurement is a linear ideal 
(non-strong) one. More precisely, let X m be the typical 
spacing between the eigenvalues of X; the measurement 
is weak when the coupling constant satisfies XX m <C A, 
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with A 2 the initial variance of the readout. One can dis- 
tinguish two regimes: the weak incoherent measurement, 
when k <§C XX m <C A, with n the coherence scale of the 
detector (in the readout basis), and the weak coherent 
measurement, when XX m <C k. The former case bears 
little interest: because of the large variance, the readout 
of the detector in each individual trial is not necessar- 
ily XX and can lie well outside the spectrum, but after 
averaging over many trials, this effect washes out, even 
if one post-selects the system [19, 20]. The latter case 
was shown to produce a large average output after post- 
selection [21, 22], and to allow a joint measurement of 
non-commuting observables with optimal noise in both 
outputs [17, 23]. It is interesting to note that a measure- 
ment of a continuous variable (X m = 0) is always a weak 
coherent measurement, and thus may show quantum co- 
herence effects when followed by a post-selection. In the 
following, no specific assumption about A are made, in 
order to have general results. For simplicity an instanta- 
neous interaction g(t) = S(t) is considered. 

In lack of any sufficient reason to believe otherwise, the 
probe and the system are assumed to be initially uncor- 
rclatcd, so that their state immediately before the inter- 
action is _R_ = ps ® f>i- We shall indicate by pi{J,J') = 
{J\pi\J') \p s {X,X') = (X\ps\X')\ the elements of the 
probe (system) density matrix in the J (X) basis and 
by = ($\piW) [ps(K,K') = {K\p s \K')\ the 

elements in the $ (K) basis. The system-probe state 
immediately after the interaction is 

Rs,f(X,X';J,J')=p s (X,X') Pi (J-XX,J'-XX'), (2) 

where the equality exp [— iJo$]| J) = | J — Jo) was used. 
After tracing out the system, the final state of the probe 
following the interaction is 

p f (J,J') = J dii{X)p s (X,X) Pi (J-XX,J'-XX), (3) 

with p{X) the spectral function describing the distribu- 
tion of eigenvalues (for a discrete spectrum, it is a com- 
bination of Dirac S-s, for a continuous one generally it 
is dp{X) = dX). The probability distribution II/ of the 
readout J is then 



n f (j) = Pf (j,j)= I d^{x) P s{x,x)Ui{j-xx), (4) 

with IT its initial distribution, and the corresponding 
characteristic function is 

Z f {4>) = [ dJe^ J n f U) = ZsiXfiZM). (5) 



Equation (5) reveals that the contribution of the detec- 
tor to the cumulants (logarithmic derivatives of Zf) are 
simply additive. In particular, if the detector is initially 
prepared in a Gaussian state, the cumulants of the output 
starting from the third and higher reflect faithfully the 



cumulants of the system, and only the second cumulant, 
which is the variance, includes a contribution from the 
detector, while the first one, the average, has a contribu- 
tion from the probe only if this is biased, introducing a 
systematic error. 

Moyal quasi- characteristic function. The results will 
be specialized to a one-dimensional system, so that X 
represents its coordinate, and P = KK its momentum. 
A useful transform of the density matrix was introduced 
in Ref. [9] : the Moyal quasi-characteristic function which 
is but the Fourier transform of W(K, X), the Wigncr 
quasi-probability function [4, 9], and is defined by 

M(x, k)=( dXdK e tkx+ixK W(K, X) 



(0) 



=JdXe^ x p{X + ^X- X -) 
= jdKe iK *p{K-^,K+\). 



From the definition Eq. (6) one realizes that M(x, 0) is 
the characteristic function for the probability Iif{K) = 
(K\p\K), and M(0, k) the characteristic function for the 
probability 11/ (X) = (X\p\X). The generalization to a 
higher dimension is straightforward. The Moyal quasi- 
characteristic function uniquely determines the density 
matrix of a system, and vice versa. For composite sys- 
tems, the marginal quasi-characteristic function of a sub- 
sytem is obtained by putting the coordinates of the re- 
maining subsystems to zero in the total function. Fur- 
thermore, the Moyal function can be expressed as the 
average of the non-Hcrmitian Weyl operator 



M(x,k) = r Yiipe ixk+ikx \ 



(7) 



where we used in the second line of Eq. (6) \X — x/2) = 
cxp[iKx/2] \X), exploited the cyclic property of the trace, 
and applied the Campbcll-Baker-Hausdorff formula. 

The usefulness of this transform shows, e.g., when con- 
sidering the measurement illustrated in the previous sec- 
tion. The joint Moyal quasi-characteristic function for 
system and probe in terms of the initial ones is simply 

M S j{x, k; = M s (x, k + X4>)Mi{<t>,j - Xx). (8) 

Results. A system interacting in rapid sequence with 
two probes, one coupling to X, the other to K, is con- 
sidered, so that the interaction term is 



Hi, 



h X x 8(t + e)$xX + X K 8(t-£)$ K K . (9) 



For e —> 0~ a measurement of K is followed by a mea- 
surement of X, and vice versa for e — > + . For e = the 
measurement is a joint measurement a la Arthurs and 
Kelly [23], that we shall not study in detail here, and 
for which we refer the reader to our previous paper [17]. 
The initial state is assumed to be R- = ps ® Pi 1 with P i 
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the density matrix of the two probes immediately before 
the first interaction. The possibility that the probes are 
initially in a correlated state is accounted for. The quasi- 
characteristic function for the system and the two probes, 
after some straightforward calculations [24] , turns out to 
be 

M s ,f(x, k- <f>,j) =M t [</>, j + A (2a s + a e A(f>)} 

xMsis + Acj)), (10) 

where for conciseness the two probes's coordinates were 
arranged in two column vectors <j> = ((/>k,4>x), j = 
(ji<,jx)i s = (x,k) represents the symplectic coordi- 
nates, we introduced the 2x2 matrix A = diag(A^, Ax), 
£ G {+, -,0}, and 







-1\ a+ + a- 

> "o = ^ • (11) 



The case e = + corresponds to a measurement of X 
followed by one of K, e = — to the opposite order, and 
e = to the joint measurement. Then, the (proper) char- 
acteristic function, i.e. the Fourier transform of the joint 
probability of observing J = {Jr, Jx) as the output, is 
obtained by tracing out the system (k = 0,x = 0), and 
by putting j = 0, so that the experimentally accessible 
quantity 



Z f {4>)= / dJ K dJ x e iJ *Il f (J K ,J x ) 



(12) 



is expressed in terms of the product of the target Ms and 
the known Mi 



Z f (<t>) =M s (A</,)M l (0,Aa £ A</>). 



(13) 



Equations (10) and (13) are the main results of this pa- 
per. Since performing a joint measurement may be more 
difficult than a sequential one, we concentrate on the 
cases £ = ±, but in principle the joint measurement of 
position and momentum allows as well the reconstruction 
of the initial state. 

If the probes are initially in a known state, then the 
unknown initial state of the system can be evinced after 
dividing the left hand side of Eq. (13) by Mi(<fi, a e ). This 



can be done as far as Mi((p, a e ) ^ 0. (A partial recon- 
struction of the state of the system may be satisfactory, 
however, if, e.g., the function M is known on a dense set 
over R 2 or everywhere but on a zero-measure set. Sec 
Ref. [25] for a detailed study.) The density matrix and 
the Wigner function are given respectively by 



W{K,X) 



dk e- lXk Z f {Vs) 
2nMi(Vs,Aa e sY 
dx e- iKx Z f (Vs) 
2?r Mi(Vs, Aa e s) ' 
dkdx e-^ Kx+x ^Zf{Vs) 



(14) 
(15) 
(16) 



(2tt) 2 M t {Vs,Aa £ s) 

with V = A- 1 = diag(A K \ A^ 1 ). 

For an ideal strong measurement, the ini- 
tial pointer density matrix would be pi(J,J') = 
5,j K ,j> K 5j x ,r x 5(J K )S(Jx), and hence Mi(<f>,j) = 
5j K ,a8j x ,o, so that the procedure would not work for 
£ = +, 4>k ^ 0, nor for £ = — , cf>x ^ 0, nor for 
£ = 0, <Pk 7^ or <f>x ^ 0. It is therefore desirable not 
to work in the strong regime, but at the same time one 
does not need to keep the measurement in the weak 
regime. For instance, the probes could be prepared in 
the mixed gaussian state 

n T n cxp{-E a [J a 2 /2Ag+ f a /2 K l}} 

with J a = (J a +J' a )/2, j a = (J a -J' a ), and k ; representing 
the coherence scale (that satisfies n a = 1/A a , where A a 
is the spread of the conjugate variable $ a , so that by the 
uncertainty principle n a < 2A a ). Then 



The ideal strong case is obtained for n a — > 0, A Q — > 0. In 
the ideal weak coherent regime, n a — » oo, A a — > oo, the 
state is Mi(<j),j) — > 6 ( j >Kt od^ >x fl (which is but the strong 
regime for when the $- variables are used as a readout), 
hence Eq. (13) can not determine Ms(\k4>k, ^x<ftx) for 
4>K,<f>x 7^ 0. Thus the procedure suggested here works 
with intermediate measurement strength. 

In particular, when the detector is prepared in the state 
given by Eq. (18), the Wigner function in terms of the 
characteristic function of Eq. (12) is 



J 



W(K, X) = 



^Z f (Vs) expJ -%{Kx + Xk) + \ 



Ak 
Xk 

Ax 
Xx 



(l + £)X x 



2 k 



x 



{l-£)X K 
2k k 



(19) 
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After substituting Eq. (12) into Eq. (19), one should be 
careful not to exchange recklessly the order of integra- 
tion, or an artificial divergence appears. Furthermore, 
for e = and n a = 2A a , A X A K = \ x \k/4, Eqs. (19) 
and (12) give a relation between W(K, X) and H/(J) 
that is the formula relating the Wigner and the Q func- 
tion, so that, as is well known [26], the joint measurement 
of position and momentum provides directly the Q func- 
tion, provided the detectors are properly prepared. By 
contrast, the general procedure proposed here for mea- 
suring the Moyal M function allows much more flexibility 
in terms of detectors preparation. 

Finally, an alternative application of Eq. (13) could 
consist in obtaining the quasi-characteristic function of 
the system for two fixed values of <pK, <Px and then vary 
the coupling strengths A_r- , Ax by keeping A^Ax fixed in 
order to reconstruct the state. 

Estimates. Let us summarize the steps needed to pro- 
ceed to the quantum state tomography with the method 
proposed here, which we may call quantum sequential 
tomography (QST), while comparing them to the anal- 
ogous steps done in quantum homodyne tomography 
(QHT). 

In QST, first, a joint probability IIj(Jx, Jx) is mea- 
sured by observing the J variables of the probe. Exper- 
imentally, one should divide the Jk , Jx plane in a suffi- 
cient number of small bins, each having an area A Jx A Jk 
determined by the precision of the measurement. The 
measurement is repeated a large number of times N, so 
that the resulting histogram approximates the true prob- 
ability within the precision of the probes. Analogously, 
in QHT, the conditional probability pr(x ( p \<p) is measured 
for different values of </> G [0, it]. The binning is in the x, 4> 
strip, and it is determined by the precision with which 
can be measured and 4> controlled. Then, in QST, the 
Moyal function Ms(x, k) is obtained by making a double 
integral, precisely a Fourier transform of the observed 
probability, which can be done efficiently fast with the 
Fast- Fourier Transform algorithm, and by dividing finally 
by the known state of the probes, which is computation- 
ally trivial. In QHT, the Wigner function is obtained by 
making an inverse Radon transform of pr(x,p\(j}), which 
consists as well in a double integral. Finally, the density 
matrix can be obtained by making a single Fourier trans- 
form of the Moyal characteristic function in QST, and of 
the Wigner function in QHT. A further Fourier trans- 
form is needed to obtain the Wigner function in QST 
and the Moyal function in QHT. Thus we can say that 
QST is more efficient of QHT for the determination of the 
Moyal function, at least as performant as QHT for the 
density matrix and less efficient for obtaining the Wigner 
function. 

Finally, we estimate the relative uncertainty in the 
Moyal function, M s (s) = Z f (Vs)M~ 1 (Vs,Aa 6 s). As 
Mi is fixed, it carries an uncertainty 5Mi(cf>,j) that 
depends on how the state of the probes was deter- 



mined. The finite sampling introduces a statistical error 
in Zf(<p), that can be estimated according to standard 
statistical analysis as SZ 2 = [1 - \Z f {4>)\ 2 ]/ N . There 
is also a numerical error introduced by the integration, 
77num- Notice how these uncertainties are present also in 
QHT. In conclusion, in QST, the relative error can be 
estimated as 

S\M s (s)\ 2 l-\Z f (Vs)\ 2 SlM^Vs.Aa.s)] 2 
\M s (s)\ 2 ~ N\Z f (Vs)\ 2 + \Mi(Vs,Aa e s)\ 2 

+ vL m (Vs). (20) 

Conclusions. A procedure for determining an unknown 
quantum state was proposed. Two measurements of con- 
jugate variables are made in quick sequence, the joint 
probability is estimated from the collected data, it is 
Fourier-transformed to give the characteristic function, 
then divided by the quasi-characteristic function of the 
probes appearing in the right hand side of Eq. (13). 
This yields the Moyal quasi-characteristic function of 
the system. The density matrix is obtained by Fourier- 
transforming the latter function. On the other hand, 
it is sufficient to determine the Moyal function only in 
a neighborhood of x — 0, k = in order to estimate 
the cumulants. Furthermore, the method proposed has 
the advantage of requiring one fixed setup and does not 
require a sharp measurement of either position or mo- 
mentum, rather it thrives over the unsharpness of the 
measurement. 

In perspective, it would be interesting to extend the re- 
sults to finite-dimensional Hilbert spaces, for which there 
is a wide interest (see, e.g., the recent Ref. [27]), espe- 
cially in the light of a recent generalization [28] of the 
concept of conjugate variables. 
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SUPPLEMENTAL MATERIAL: 
DERIVATION OF THE MAIN RESULT. 



The case of a measurement of momentum followed by one of position is considered for dcfmitcncss. The coupling 
constants are temporarily absorbed in the rescaled variables $ a — > A a $ a , J a — > J a /A a , a = K,X. The joint density 
matrix after the interaction is 

R SJ = e i*xX e i* K K ft_ e -i& K K e -i*xX ^ (21) 

with i?_ = ps ® Pi the initial joint density matrix. The most convenient representation is in terms of |$) for the 
probes and |X)for the system 

R SJ (<f, X, X') = e^ x Ps (X + $ K ,X' + & K ) Pi ($, (22) 

where exp[— iaK] \X) = \X + a) was used. Multiplication by exp[ikX] and by exp[ij • <£>], with X = (X + X')/2 and 
<§ = ($ + $')/2, followed by integration over X and $ yields 

M SJ (x, k; = J d<bdX e *M*+**+<*+**>*] ft ($ + ± t $ - ^)p s (X + * K + ^!L,X + # K - £±^) 

d i e .Mx+^-(Hfc)* K ] ft ($ + t } $ _ t) Ms {x + 4> K , k + 4> x ) 

= Mi(tj>,j + l)M s {x + 4> K , k + (f> x ) (23) 

with 7 = (—k — 4>x, x). The identity <&xX — ^' X X' = X<fix + &xx was used, and the change of variable Y = X + <&k 
was made. 
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